PLOS Biology

Integrating multisensory and behavioral information is essential for sensory perception. In the auditory system, multisensory and behavioral influences emerge early in subcortical structures. Descending projections from non-auditory cortical areas are well positioned to convey such signals, yet how they shape subcortical auditory processing remains poorly understood. Here, we investigated corticocollicular projections from the primary somatosensory (S1) and motor (M1) cortices to the inferior colliculus (IC), a principal integration center in the auditory midbrain. We found that trunk- and limb-related regions of S1 and M1 form prominent monosynaptic projections to the IC, and that optogenetic activation of these projections robustly drives IC activity. Notably, a substantial population of cortical-responsive neurons did not respond to sound. In sound-responsive neurons, concurrent cortical stimulation enhanced sound-evoked responses, whereas cortical activation preceding sound onset suppressed them. Furthermore, both cortical-responsive IC neurons and deep-layer S1 and M1 neurons exhibited locomotion-related modulation and anticipatory activity prior to movement onset, suggesting that these descending pathways convey movement-related signals to the IC. Together, our findings identify a descending sensorimotor circuit that integrates body- and movement-related information with auditory processing in the auditory midbrain.

Background Using a population-based cohort from the Calgary Health Zone (CHZ), Canada, we integrated longitudinal antimicrobial susceptibility and prescribing data with the whole genome sequences of five major pathogens. We aimed to assess how antimicrobial resistance (AMR) responds to prescribing changes and determine which bacterial strains shape these dynamics. Methods We analysed antibiotic prescribing rates, clinical and genomic data from 7,271 Staphylococcus aureus, 1,609 Enterococcus faecalis, 801 Enterococcus faecium, 11,363 Escherichia coli, and 2,319 Klebsiella pneumoniae isolates, associated with bacteraemia episodes in the CHZ between 2006-2022. Genomic clusters (referred to as strains) were identified using StrainGST and assigned to known sequence types (STs) or clonal complexes (CCs). Strain-level incidence, stratified by community-onset (isolates collected [&le;]48h after admission) and hospital-onset (>48h after admission), AMR phenotypes, and prescribing rates were modelled using negative-binomial and binomial regression. Temporal trends were quantified using average annual percentage change (AAPC). Findings Between 2010-2022, fluoroquinolone prescribing declined in both community (AAPC=-6.8% [95% CI -8.1, -5.4]; p<0.0001) and hospital settings (AAPC=-5.1% [-6.5, -3.7]; p<0.0001). This was accompanied by a significant reduction in fluoroquinolone resistance among Gram-positive species. Specifically, S aureus bacteraemia resistant to clinically important antibiotics, cloxacillin, ciprofloxacin, erythromycin, and clindamycin, declined from 2006 to 2022, mostly in hospital-onset cases (AAPC=-16.0%, [-19.3%, -12.7%], p<0.0001). In E coli, ceftriaxone and ciprofloxacin resistance were clustered in ST131 and the emerging ST1193; the latter increased steadily, particularly in community-onset cases (AAPC=17.7%, [0.0%, 30.0%], p<0.0001). CTX-M-27-producing E coli ST131 strains increased (AAPC=23.8%, [17.4%, 30.5%], p<0.0001) between 20082022, while CTX-M-14-producing E coli ST131 declined (AAPC=-15.9%, [-21.3%, -10.2%], p<0.0001) between 2013-2022. These trends were paralleled by an increase in community cephalosporin prescribing (AAPC=7.3%, [4.2%, 10.5%], p<0.0001) between 2010-2022. For K pneumoniae, hypervirulent ST23 was most common (N=88) with an increasing trend in incidence (AAPC=3.0%, [-2.8%, 9.2%]) between 2006-2019. Conclusions The contrasting resistance trends between Gram-positive and Gram-negative species underscore the complexity of AMR control efforts. Effective strategies will require stewardship efforts targeting multiple drug classes, genomic surveillance for emerging resistant strains, and interventions extending beyond hospital settings.

Reputation plays a major role in supporting cooperation among unrelated individuals through indirect reciprocity. By helping others, individuals build a good personal reputation and receive greater benefits from future partners. Most models of indirect reciprocity assume that a person's reputation reflects only their own behaviour. Yet in many societies, people are also judged by their family's reputation. How family reputation affects the evolution of cooperation, and whether reliance on it can itself evolve, remain unclear. Here we show that reputation inheritance expands the conditions under which indirect reciprocity favours cooperation, increasing helping and favouring greater reciprocity. Greater reciprocity in turn favours stronger reliance on inherited reputation, creating a positive feedback that stabilises cooperation, especially when interactions are infrequent or personal behaviour is difficult to observe. This feedback arises because cooperation generates future benefits both for the individual, through their personal reputation, and for their descendants, through inherited reputation. Reputation inheritance thereby provides a route via which kin selection and reciprocity, often treated as alternative explanations for cooperation, can reinforce one another. Our model helps explain why family-based reputation occurs across diverse human societies and provides an evolutionary framework for studying phenomena organised around family standing, including kin-based institutions, feuds between families and honour-based violence within them.

The evolution of new traits can drive species diversification by facilitating the use of new resources, but environmental change may turn these same adaptations into liabilities.Trait loss is also often associated with the origin of new ecologies, but how losses modulate diversification remains unclear. The swim bladder allows ray-finned fishes to regulate their buoyancy and exploit ecosystems throughout the water column, yet this organ has been lost many times among species-rich lineages. Here, we show that timing and ecological context control the macroevolutionary effects of swim bladder loss. Many lineages of fishes lost the swim bladder over the last 66 million years as they specialized for benthic habitats where buoyancy regulation is unnecessary. Swim bladder loss enabled the descendants of these benthic fishes to diversify in the deep sea where extreme pressure makes its inflation untenable, and in the frigid, oxygen-saturated Southern Ocean, where loss of the oxygen delivery mechanisms required for swim bladder inflation carries little physiological cost. Yet, we detect a selective filter associated with swim bladder loss during extreme global warming 56 to 50 million years ago, when its absence limited the capacity of fishes to escape ecological disruptions on the ocean floor. These contrasting patterns explain how the loss of a complex trait promoted major ecological transitions without increasing overall diversification through deep time. As human activity drives rapid global warming, the evolutionary legacies of swim bladder loss may again shape the fate of marine fish diversity.

Animals rely on color to navigate complex environments, yet how eyeless organisms use chromatic information to guide food choice remains poorly understood. Here, we show that Caenorhabditis elegans exhibits robust color dependent foraging driven by microbial chromophores, preferentially consuming red while avoiding blue chromoprotein expressing bacteria across bacterial backgrounds and wild isolates. This discrimination persists in darkness and independently of photoreceptor, revealing a mechanism beyond canonical photoreception. Purified chromoproteins and bacterial metabolite fractions independently reproduce preference, demonstrating complementary chromatic and post ingestive metabolic cues. Mechanistically, blue chromoproteins generate singlet oxygen, producing oxidative stress and remodeling bacterial tryptophan and pterin metabolism, whereas red food promotes serotonin production and feeding-associated neuropeptide signaling. Disrupting serotonin biosynthesis or neuropeptide processing abolishes color preference. Together, our findings reveal a previously unrecognized, novel sensory strategy in which wavelength-selective pigment photochemistry transforms microbial color into metabolic information that is integrated through gut brain neuroendocrine signaling to guide foraging behavior in an eyeless animal.

Microorganisms across all domains of life engage in molecular conflict, deploying toxins to inhibit competitors or respond to biological threats. Among these, ribonuclease toxins are particularly widespread and diverse. A substantial fraction is associated with the BECR fold, a compact /{beta} architecture that supports RNase activity despite extensive divergence. Although several canonical members are well characterized, many BECR-fold proteins remain difficult to identify because of low sequence similarity, variation in catalytic residues, and structural elaborations that obscure evolutionary relationships. The growing availability of high-confidence protein structure predictions provides an opportunity to reassess this deeply divergent protein landscape. Here, we integrate iterative profile-HMM searches, profile-similarity networks, structural analyses, active-site mapping, and genomic context to examine BECR proteins across the tree of life. Our analysis resolves an expanded BECR-fold landscape comprising canonical BECR and BECR-like superfamilies, refines the organization of canonical BECR proteins and identifies previously unrecognized families. We further validate BECR-Tox2 as a toxin neutralized by a cognate immunity protein and show that its homologs occur in both Menshen-like anti-phage systems and polymorphic toxin loci. Together, these findings expand and clarify the BECR-fold landscape and provide a framework for identifying and interpreting highly divergent proteins of this fold.

Phage D29 infects Mycobacterium smegmatis mc2 155 and has a non-canonical lysis cassette that encodes two endolysin proteins (Lysin A and Lysin B) and a single two transmembrane domain (TMD) protein, LysA2a similar to F1 cluster phage LysF1a. A 1TMD LysF1b homolog, LysA2b, is encoded by a gene found downstream of the tape measure. Exogenous expression of both LysA2 proteins in tandem is a cytotoxic to M. smegmatis. Deletion of lysA2a produces phages that are lysis competent with a 10-minute triggering delay and 30% plaque size reduction. Deletion of lysA2b results in severe lysis defects manifest by 70% reduced plaque size, delayed lysis timing and reduced burst size. Deletion of both lysA2 genes results in phages that are viable and show lysis phenotypes like the lysF1b deletion. Genetic complementation of lysA2b deleted phage with the lysF1b gene fully complements the lysis phenotypes but alters the triggering time to that of an F1 cluster phage. Energy poisons trigger lysis prematurely in all phages with lysA2 gene deletions. Lysis recovery mutants (LRM) isolated from phages lacking the lysA2b genes generate wild type plaque size and have point mutations that map to TMD1 or the C-terminal region of the lysA2a gene. LRMs isolated from phages lacking both lysA2 genes show premature lysis and have mutations that all map to residue C31 of a novel lipoprotein (gene 64). Deletion of gene 64 does not change wild type D29 lysis phenotypes or rescue the lysis defects of any of the lysA2 mutants. A fitness/competition assay shows that loss of the lysA2 genes imposes a substantial competitive fitness cost. These finding support a lysis regulatory network model where the 2TMD protein is maintained in an inactive state until activated by its cognate 1TMD lysis regulator and the lipoprotein has accessory function that may enhance lysis efficiency.

A dental AI system and a dentist reading the same radiographs form a paired comparison. Published comparative studies often report the two arms separately against a reference standard, leaving the joint pattern of correctness between them unavailable for secondary paired inference. We show what that omission costs. The accuracy difference remains exactly identified; its sampling variance does not, so the report contains the estimate and not its uncertainty. On a study of 282 units, two published accuracies are consistent with 38 distinct joint tables whose confidence intervals differ in width by a factor of 2.5. The consequence is a three-zone decision map rather than a single threshold: differences at or below 1.06 points are non-significant under every compatible table, differences at or above 6.03 points are significant under every compatible table, and in between the published numbers cannot decide. We then show the omission is repairable at negligible cost. One additional integer, the number of units both arms classify correctly, identifies the joint table exactly and restores standard paired inference. For a panel of readers the pairwise dependences must arise from one joint distribution, a constraint that binds once three readers are present; publishing each reader's joint-correct count against a single reference reader cannot widen and may tighten every pairwise bound, and in a 7-arm experiment reduced them by a median of 37% even for pairs excluding that reference. Where the integer was never published we give DentalPair-Cert, an interval with finite-sample coverage uniformly over every admissible within-unit AI-dentist dependence under the independent-sampling-unit model, certified in both the nuisance maximization and the inversion. Across 4,200,000 simulated comparisons an independence analysis falls to 74.5% coverage with 12.2% type-I error; in a purposive sample of 9 recent comparative studies, 1 reported a paired test on discordant units.

Machine-learning models of the human microbiome are trained mostly on stool samples from single cohorts, limiting cross-body-site representation and cross-study generalization. Progress is constrained less by algorithms than by the absence of a harmonized multi-body-site corpus carrying the technical metadata needed to model, rather than ignore, batch structure. Here we release Corpusome, a harmonized two-tier cross-body-site human microbiome corpus for representation learning: a harmonized corpus of 187,546 human microbiome samples integrating standardized profiles from curatedMetagenomicData, the American Gut Project, and the EBI MGnify platform. Corpusome follows a two-tier design preserving both functional depth and cross-body-site breadth: a shotgun tier (22,588 samples, 93 studies) with species- and pathway-level profiles, and a 16S tier (164,958 samples, from a full pull of 708 MGnify studies) with genus-level profiles extending coverage to oral, skin, respiratory, and urogenital sites. It spans six body sites and two modalities, with harmonized metadata for batch-aware modelling. Body-site signal exceeds technical/source variance in the 16S tier by approximately 2.4-fold.

Next-generation sequencing (NGS) of the prokaryotic 16S rRNA gene revolutionized gut microbiome research two decades ago. However, short read lengths remain an inherent limitation of platforms such as the widely used Illumina platforms (2 x 150-300 bp). Recent advances in Oxford Nanopore Technologies (ONT) flow cell chemistry (R10.4.1) have substantially improved sequencing accuracy. Combined with a custom multiple-primer strategy that comprehensively targets 16S rRNA gene variants to generate near-full-length amplicons, this approach enables read-by-read taxonomic classification, a feature not feasible with short-read sequencing platforms. Although our multiple-primer strategy could enable parallel sequencing of more than 18,000 samples (192 x 96), current flow cell capacity offers sufficient sequencing depth for approximately 1,000-1,500 samples. To validate the scalability and our per-read classification pipeline, we show that more than a thousand human fecal microbiome samples spiked with two bacterial strains (Imtechella halotolerans and Allobacillus halotolerans), not otherwise present in human fecal samples, can be successfully sequenced on a single flow cell, achieving a per-molecule error rate sufficient for direct per-read classification and at an adequate read depth for downstream analysis. This level of scalability significantly reduces per-sample costs, making the approach more accessible to a broader research community. To embrace these advancements, we have developed RubyRed, a pipeline that processes raw sequencing data and assigns taxonomic classifications on a per-read basis. Using spike-in references (I. halotolerans and A. halotolerans), we demonstrate high mean single-read sequencing accuracy (99% and 98.9%, respectively), with the majority of reads exceeding the canonical threshold required for species-level taxonomic classification based on the 16S rRNA gene.

Motivation: The Interactive Tree of Life (iTOL) is widely used to display and annotate phylogenetic trees, but managing its format-sensitive annotation files impede reproducible high-throughput analyses. Among the maintained Python packages and versions evaluated, none combined template generation, taxonomic monophyly assessment and iTOL batch operations. Results: PyiTOL validates inputs, generates 31 iTOL template schemas (22 accepted by the live batch uploader), performs LCA-based monophyly classification with nested-monophyly detection, sampling-completeness states and polyphyletic subgroup decomposition, plus API upload and session replay. On a topology-constructed benchmark, all calls matched prespecified labels for 4,389 groups; on a 700-genome tree, binary mono/non-mono calls agreed with ETE4 for 409 genera; 17,294 GTDB R232 genera were processed in about 17 s. Availability and Implementation: PyiTOL 1.0.3 (Python [&ge;]3.10; Linux, macOS and Windows) is MIT-licensed at https://github.com/ZengZichao/PyiTOL and archived with test data at Zenodo (https://doi.org/10.5281/zenodo.22106806).

Pneumococcal vaccination of older adults is primarily guided by age and clinical eligibility, despite substantial variation in individual risk of severe pneumonia. Here, we used longitudinal electronic health records from 787,538 adults aged [&ge;]65 years to evaluate the real-world effectiveness of the 20-valent pneumococcal conjugate vaccine (PCV20) and quantify clinical benefit according to baseline risk of pneumonia hospitalization. We developed and validated a machine-learning model using pre-PCV20 data to estimate individual 12-month hospitalization risk and integrated these predictions into a propensity score matching framework. Overall vaccine effectiveness against pneumonia hospitalization was 16.5% (95% CI, 10.6-22.1), but this population-level estimate masked substantial heterogeneity in clinical benefit. The 60% at lowest predicted risk, characterized by younger age and fewer pulmonary and other chronic conditions, showed no measurable reduction in hospitalization (VE, 3.1%; 95% CI, -14.4 to 18.0) and had an estimated 1-year number needed to vaccinate (NNV) of 7,423, compared with 184 and 115 in the intermediate- and high-risk groups, respectively. These findings suggest that incorporating baseline risk into adult pneumococcal vaccination strategies could enable more targeted and potentially better-timed vaccination.

Vaccine effectiveness (VE) estimates can vary widely between years and populations, even for the same vaccine. Estimated VE is known to be sensitive to susceptible depletion and differences in pre-vaccination infection risk between vaccinated and unvaccinated populations. However, how variation in pre-vaccination risk within and between the two groups affects VE estimates over time remains unclear. This uncertainty is especially important given negative VE estimates. We investigated the difference between estimated VE and true vaccine protection considering continuous distributions of pre-vaccination infection risk under three scenarios. When the vaccinated and unvaccinated populations differ in their mean risk, estimated VE can be higher or lower than true vaccine protection. Similar patterns arise when both populations share identical means but different risk distributions. Finally, if infection-derived immunity lasts longer than vaccine protection, annual VE estimates can vary by tens of percentage points between years despite constant true vaccine protection. These theoretical results underscore that VE studies estimate contrasting risk between vaccinated and unvaccinated individuals in a particular time and place, and VE estimates can vary counterintuitively between years and populations even with constant vaccine-induced protection. Explaining variability in estimated VE thus requires a more complete understanding of populations' distributions of infection risk.

Autosomal dominant gain-of-function mutations in platelet-derived growth factor receptor beta (PDGFRb) cause overgrowth of the skeleton and other connective tissue in Kosaki overgrowth syndrome. However, the target cell type and signaling pathways underlying PDGFRb-driven overgrowth are unknown. Normal postnatal growth is controlled by pituitary-secreted growth hormone (GH), which activates the STAT5 transcriptional factor to upregulate insulin-like growth factor 1 (IGF1). To investigate the role of the GH-STAT5-IGF1 pathway in PDGFRb-related overgrowth, we generated mice with a PDGFRb gain-of-function mutation in skeletal and fibroblast lineages, which resulted in STAT5 activation and gigantism. Conditional deletion of Stat5ab in connective tissue lineages rescued skeletal overgrowth and keloid-like fibrosis in the skin. Conditional deletion of GH receptor (Ghr) did not rescue overgrowth, indicating the physiological activator of STAT5 is not required for overgrowth. However, deletion of Igf1, the STAT5 target gene, and its receptor, Igf1r, in connective tissue, rescued the overgrowth phenotype. These findings demonstrate a GHR-independent STAT5-IGF1 signaling pathway in mutant connective tissue cells, which mediates PDGFRb-driven overgrowth in mice and potentially in humans with similar PDGFRB mutations.

Post-acute infection syndromes (PAIS) follow viral syndromes including post-acute sequelae of COVID19 (PASC) which complicates 10-25% of SARS-CoV-2 infections. These syndromes lack precise explanatory mechanisms. We studied 173 human saliva proteomes during respiratory viral syndromes, seeking associations between 44 clinically-relevant protein expression patterns and subsequent sequelae counts. Exploratory models adjusted by extensive clinical annotations found interactions between 23 acutely-responsive proteins and SARS-CoV-2 infection that inversely predicted subsequent neurocognitive sequelae. An overlapping 19 acutely-responsive proteins during any acute respiratory viral syndrome inversely predicted general fatigue-related sequelae. Altogether, 29 proteins, derived from interferon stimulated genes (ISG), were uniformly beneficial, including 13 predictive of both neurocognitive and general sequelae. The proteins suggested both shared early pathobiology and virus-specific protective responses that shaped resolution of acute disease and different PAIS. Acutely elevated protective ISG proteins associated with reduced post-viral symptoms identify investigational starting points for novel mechanisms, diagnostics and therapeutics for PASC and PAIS.

Forecasting influenza hospitalizations informs public health preparedness, yet questions remain about which types of forecasts best guide action. We evaluate categorical trend forecasts, which communicate probabilities of upcoming increases or decreases in epidemic trajectories, submitted to CDC's FluSight Forecasting Challenge between Fall-2024 and Spring-2026. Teams submitted probability distributions over five categories describing direction and magnitude of week-over-week changes in laboratory-confirmed influenza hospital admissions. We assessed performance using Ranked Probability Skill Score, Brier Skill Score, and measures of forecast-observation agreement. Most models outperformed an equal-probability baseline; the FluSight ensemble ranked among the top three in the 2024-25 and 2025-26 seasons. Forecasts were most accurate during stable periods and least during periods of rapid change, with most models underestimating observed trends. Conclusions were robust to choice of scoring metric and reference model. These results support categorical trend ensembles as an approach to communicating infectious disease forecasts that may inform public health decision-making.

Background Understanding the mechanisms underlying brain aging and age-related pathological changes is essential for advancing brain health research. Our group previously developed a mechanistic mathematical model of healthy brain, Chamberland et al. (2024) that integrates key biological processes involved in normal aging, from which Alzheimer's disease (AD) related changes may emerge naturally. Objectives To characterize and validate this brain model by evaluating its sensitivity, calibrating its parameters, and assessing generalizability in independent populations. Methods The model represents the evolution of key biological processes associated with brain aging, including amyloid beta (A{beta}), tau pathologies, neuroinflammation, and neuronal death. After identifying the 30 most influential parameters, we calibrated the model using cognitively normal (CN) participants from the AD Neuroimaging Initiative (ADNI) database (n = 211) by minimizing a loss function composed of three outcomes (AB) plaques, tau tangles, and neuronal density). The calibrated model was then applied to the UK Biobank cohort (n = 35,899) of normal controls (aged 44-82 years). The effects of sex and APOE were evaluated using stratified simulations. Results Parameter calibration significantly reduced the prediction errors for A{beta} and tau. Neuronal density predictions showed strong agreement in the UK Biobank cohort. The variance decomposition identified APOE status as a major contributor to variability in A{beta}. Conclusion Our validated brain health model links mechanistic pathways with population data and reproduces neuronal density patterns in an independent cohort. These findings support its use as a framework for studying brain aging and investigating how Alzheimer's disease related pathological changes may emerge with aging.

Increased activity of leucine-rich repeat kinase 2 (LRRK2) confers Parkinson's disease risk. LRRK2 dynamically localizes to lysosomal membranes in response to various stresses, yet the mechanisms by which distinct lysosomal perturbations are communicated to LRRK2 remain unclear. Here, we show that inhibition of the lysosomal lipid kinase PIKfyve promotes LRRK2 recruitment and signaling through a pathway that requires the lysosomal chloride/proton antiporter ClC-7. ClC-7 in turn controls the accumulation of multiple Rab GTPases on lysosomes. LRRK2 signaling under these conditions requires its established Rab-binding surfaces, with Rab12 contributing significantly to this response. This pathway operates independently of CASM. In contrast, lysosomal stresses that induce CASM require both Rab-binding sites on LRRK2 and GABARAP for robust LRRK2 signaling. These findings identify ClC-7-dependent lysosomal remodeling and Rab accumulation as key features linking PIKfyve inhibition to LRRK2 signaling and reveal that distinct lysosomal stresses engage different combinations of Rab and GABARAP inputs to activate LRRK2.

Mood fluctuations in major depressive disorder are difficult to anticipate. The biological neural rhythms that organize mood dynamics over days to weeks remain unknown. In individuals implanted with a chronic neural sensing and stimulation device for treatment-resistant depression, we collected years-long intracranial neural recordings alongside daily mood ratings. Both mood and limbic neural activity fluctuated cyclically with multiday (multidien) periodicities of 2-34 days. An individual's daily phase position within mood cycles tracked depression severity, distinguishing whether symptoms were rising, peaking, or resolving. Neural rhythms led mood cycles and forecast an individual's mood trajectory up to 30 days in advance, outperforming models based on raw neural activity. Electrical stimulation reshaped these rhythms, shifting individuals away from the peak-depression phase of their multidien cycle. Our results identify multidien rhythms as an organizing principle of mood in depression and a forecastable, modifiable target for chronotherapeutic neuromodulation.

Wildfire smoke is an escalating global public health threat exposing millions of people, including children, to hazardous air pollution each year. Although wildfire smoke toxicants have been linked to a range of adverse health outcomes, including immune dysregulation, the long-term consequences of real-world pediatric wildfire smoke exposure on health and development remain largely unknown. To investigate the persistent effects of early-life exposure on immune health, here we leveraged a cohort of rhesus macaques that experienced nine consecutive days of hazardous wildfire smoke exposure in infancy during the 2020 Oregon Labor Day wildfires. By integrating ex vivo immune stimulations, multiplex cytokine profiling, single-cell transcriptomics, and genome-wide DNA methylation profiling, we identified persistent immunological consequences across molecular and functional levels. We found that a single severe postnatal exposure, in the first three months of life, was associated with persistent change in the innate immune response, including reduced pro-inflammatory cytokine response to a bacterial endotoxin, with subtle but consistent transcriptional changes in myeloid cells, particularly among males. Wildfire smoke exposure was also associated with changes in proportion of B and T/NK cells, and within the T/NK cell compartment, exposed animals exhibited an expansion of cytotoxic cells. Consistent with this, CD8+ T cells displayed extensive transcriptional remodeling and shifted toward more differentiated effector states, with the greatest differentiation observed in animals exposed at the youngest ages. Genome-wide DNA methylation profiling identified smoke-associated methylation changes consistent with acceleration of epigenetic aging, as well as persistent epigenetic alterations impacting genes involved in oxidative stress responses, innate immunity, T cell differentiation, and hematopoiesis. These findings demonstrate that a single severe wildfire smoke exposure during a critical developmental window is associated with extensive immune and epigenetic remodeling that persist years after exposure, providing new insight into the long-term biological consequences of early-life wildfire smoke exposure.